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S ' Abstract 

In this paper we study the orbits of massive bodies moving in the 
spacetime generated by a spherically symmetric and non-rotating dis- 
tribution of mass. More specifically, our treatment discusses the more 
accurate calculation of the precession of pericenter due to general- 
relativistic effects. Our result is accurate up to terms of second order, 
y—i , while the precession met in the bibliography is accurate only up to 

first-order terms. 

in : 

1 The Schwarzschild spacetime 

\o '. 

The first solution of Einstein's field equations was published by Karl Schwarz- 
"q \ schild in 1916. It is a solution that informs us about the spacetime generated 

in the exterior of a spherically symmetric and non-rotating distribution of 
mass. As is well known, that distribution of mass provides the region that 
surrounds it with a static, spherically symmetric spacetime. That kind of 
• £h ! spacetime is mathematically denoted by the line element (cf. [Q) 

ds 2 = g tt (r)dt 2 + g rr (r)dr 2 + r 2 d0 2 + r 2 sin 2 9 dcj) 2 

= _e2* df a + e 2A dr 2 + ^2 + r 2 ^2 g # 2 

in spherical polar coordinates (t,r,6,<p). In the second step above we intro- 
duced the functions $ = $(r) and A = A(r) in place of the two unknowns 
gu{r) and g rr (r) respectively. That replacement was possible since g u < 
and g rr > anywhere in spacetime. 1 Of course, we have to impose a couple 
of boundary conditions to the aforementioned line element, i.e. 



lim $(r) = lim A(r) = 



r^oo 
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1 Bear in mind, however, that this allegation breaks down in the case of black holes, 
where we should reconsider our system of coordinates. 
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for we demand, as is physically reasonable, that spacetime be flat far away 
from the distribution of mass. 

The previous results give us the ability to calculate the components G^ v 
(/i, v = 0, 1, 2, 3) of Einstein's tensor. The missing step, now, in order to 
calculate the unknown functions $(r) and A(r), is to find out the components 
T^ v of the stress-energy tensor, plug them into Einstein's field equations, 

G*» = 8ttT^ 

(note that here and hereafter we use geometrized units unless otherwise men- 
tioned) and solve the resulting differential equations for $(r) and A(r). The 
aforementioned tedious calculations lead us to the expressions 2 

e 2 * = e" 2A = 1 - — 
r 

where M is the total mass of the distribution of mass. Therefore, we are 
finally in position to write down the line element in its final form: 

r/.s- = ( 1 - —) dt 2 +(l- —) 1 dr 2 + r 2 d6 2 + r 2 sin 2 6 deb 2 




2 Precession of Pericenter 

In order to calculate the precession of the pericenter of the elliptic orbit of a 
massive body revolving an attractive center, we begin from the equation of 
motion 

where we suppose that the elliptic orbit is described by the radial distance, r, 
and the azimuthal angle, 0, r is the proper time, E and L are the energy and 
momentum per unit mass respectively and M is the total mass of the attrac- 
tive center. The shape of the effective potential V 2 {r) — (l — ^) ^1 + 

for a typical value of the parameter L 2 is shown in Fig. 1. Our study will 
be limited to energies E 2 « 1, for we want the body to be in a stable orbit 
around the attractive center. As an example we can see the points A and B 
in Fig. 1, among which the body can move in a nearly elliptic orbit. 

In order to write equation in the more convenient form dr /deb = f(r) 
we use the fact that 

^ = V 4> = E± = g<M>P± = L~ L 

dr m m r 2 



2 For further details see £Q 
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Therefore, we get 

dr 



2 _E 2 -(l- 2 f) (l + g) 



L 2 



We now proceed with the transformation u — —, where r is a constant with 
length dimensions. That transformation gives us a differential equation with 
no dimensions: 



du E 2 -l _ 2Mr „ 2M _ 
— = \ — ~ ri H = — u - u l H n d 



It is clear that we can separate the variables in the last differential equation, 
thus taking 

du 



E 2 ~l r 2 I 2Mr _ 2 i 2M 3 



Let, now, a, 6 and c be the roots of the polynomial 

n . . E 2 -l 2 2Mr 2 2M 3 
p( w ) = = r l -| = — --y — u -\ u 

L 2 L 2 r 

Then, P(u) can be written as 

P{u) = a(u — a)(b — u)(c — u) 
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where a = 2M/r . Because of the equality of the coefficients of the same- 
order terms for the two forms of P(u) we get the equations 



abc = ^rl 2 

2ML 2 

2 

ab + c(a + b) = ^ (3) 

and 

We suppose that within the range of values of energy for which we study 
our problem, P(u) has the three real roots, a, b and c for which a < b ^ c. 
Therefore, for a < u < b we can state that the relation c ^> u holds. That 
statement enables us to Taylor-expand the term 1 / y/c — u, 

1 1 f 1 u 



and get 



y/c — U y/c V 2c 

—^^^^^^=^^^^^^= = dtp =>- 

y/a(u — a) (b — u) (c — u) 

du 1 udu 



Vac yj(u — a) (b - u) 2cy/ac ^{u - a)(b - u) 

We integrate the last relation for u from a to b, so varies from zero to 
0finai/2. The calculations yield 



du 1 f udu 



^final /2 



olc J y/( u - a )(b- u) 2cy/ac J ^{u - a)(b - u 



1 1 a + b fina i 

= 7T H 7= — - — 7T = — — (5) 



y/ac 2cyfac 2 



\ " 6 b 



for J , = arctan 



7r and f ; = = ^7T. 



a yj (u-a)(b-u) \ V (u-a)(b-u) 

Let, now, d + 6 = e < c. From equation (jlj) we get c = ^ — e and so 

1 Me 3 M 2 e 2 
1 + + 



ac r 2 r. 



2 



4 



Therefore, equation (J5J) becomes 



-'final 



2vr + 2tt : 



Me 

To 



3vr- 



M 2 e 



2 r 2 



H 7T 1 + 

2c 1 



Me 3 M 2 e 2 



r 



+ 



2 r 2 



where we neglect terms of order higher than (Me/ro) 2 . Apparently, the 
resulting precession is 

A0 = fina l - 27T 

Me 3M 2 e 2 ' 



Me M 2 e 2 e , 

A0 PS 27T h 37T 5 1 7T | 1 + 

r 2c 



r 



2 r 2 



Me 

~ 2tt h 3tt- 



M 2 e 



2^2 



+ 



r 



2 (^- e ) 



Me 3 M 2 e 2 

7T | 1 + + 2— 

r 2 



Me 



M 2 e 



2^2 



Me 



A0 ~ 2lT h 37T 5— + 7T 



r 



r 



2Me 4M 2 e 2 

1 + + s— 

r 



Me 3 M 2 e 2 
r 2 



Me M 2 e 2 

A0 PS 37T h 67T 5— 

r 



2Me\ 



(6) 



-1 



keeping, 



r 

where we Taylor-expanded the term (t^j — e) 1 = (l — 

one more time, terms of order not higher than (Me/ro) 2 . 

In order to calculate e we will use equations (J2J) and (0J), from which we 

get 



7-q 
2M 



e e + 



£ 2 



1 



2ML 2 °^ 



2A1 



V- 



Again, if we use the Taylor expansion of the term — e) 1 which appears 
in the last equation, we take 



4M 2 (l-£ 2 ) 

E 2 



2M 



E' 



2Mr ) e + ^-E 

L 2 







The straightforward solution of the last equation gives us the values 
i-£ 2 < 

2M ~ r 

e 



rp _i 

2M ~ LP- 



2Mr 



_ 4M 2 (1-E 2 ) 
L 2 



± 



_ro_ 

2A1 



L 2 



1 _ AM 2 (1~E 2 ) 
L 2 
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From the two possible values of e (positive and negative sign) we choose the 
one with the negative sign. That choice is based on the fact that we need e to 
be a small compared to (cf. equation (jlj) and the discussion that follows 
it). The negative sign validates the requirement e<c~ something that 
is not accomplished via the use of the positive sign. Therefore, that choice 
and the fact that for E 2 » 1 it is ( E 2 + i) (1 - E 2 ) « |(1 - E 2 ), enable us 
to write that 



rp , l-£ 2 
2M L' 2 



ro 

2Af * 



ft 16Af 2 I 24Af 2 
L 2 ~^ Z 2 



^(l-£ 2 



, + 6gi(l 



•£ 2 ) 



1 _ 4M 2 (1-E 2 ) 
L 2 

If we Taylor-expand the square root above and keep terms of order not higher 

m£ 

L 4 



than y—, we get 



1 + l^Ml4M 2 
r T l 2 



1 6M 4 
L 4 



4M 



4A/f 2 (l-g 2 ) 
L 2 



where we nee; lect the term ^-(1 - E 2 ) 2 which appears during the calcula- 



tions since E 2 « 1. Finally, we Taylor-expand the term 
and we end up with 



4M 2 (1-E 2 ) 
L 2 



2ME 2 8M 3 , ~ 



(£ 2 - 3) (e 2 - i) 



r 



L 2 L 4 

If we replace the acquired value of e to equation (JHJ) we get the final result 



(7) 



A0 



6vrM 2 
L 2 



UM 2 \ 



6M 2 



in terms of the total mass of the attractive center and the energy and mo- 
mentum per unit mass. 

Moreover, we can express the precession found (equation (jHJ)) as a function 
of the eccentricity, e, and the semi-major axis, a, of the elliptic orbit of the 
body. If we assume that the mass of the attractive center is much larger than 
the mass of the orbiting body, then 



Ma(l - e 2 ) 



and 



E 2 = 1 - M/2a 



Therefore the precession in terms of e and a is given by the formula 



A0 



6vrM 
a(l-e 2 



2a V 1 - e 2 



(9) 
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where we neglect terms of order higher than (M/a) 2 . 

Equation can be used to calculate the precession of the perihelion 
of planet Mercury's elliptic orbit around the Sun. The semi-major axis of 
Mercury's orbit is a = 5.791016 x 10 10 m, its eccentricity is e = 0.205615, 
while the mass of the Sun is M = 1.9892 x 10 30 kg. Of course, in order 
to find some real numerical values we should write down equation (^J) in 
SI units. Equation (JUJ) as we see it is valid as long as we remember that 
we use geometrized units, where the speed of light, c, and the constant of 
gravitational attraction, G, are taken to be equal to one: c = G = 1. The 
form of equation © in the SI system of units is 

2ac 2 V 1 -e 2 y_ 

where G = 6.672599 x 10- n N-m 2 /kg 2 and c = 299792458 m/s. So if we take 
into account the fact that Mercury spins around the Sun 415 times within a 
century, then we get the precession 

A^Mercury = 42.964926" per century 

Note that in common bibliography the term 1 — (l — y^-) — whose nu- 
merical value is 1.000004813 in the case of Mercury — is just equal to unity. 
That exlains the fact that the precession we meet in the bibliography is 
A^Mercmy = 42.964720" per century. It is clear that the precession derived 
by this new approach for the case of planet Mercury does not differ much 
from the precession derived with the usual way pQ. However, if we happen 
to study the precession of pericenter for an orbit near a strong gravitational 
field (e.g. in the vicinity of a black hole), then the extra term is expected to 
produce additional precession that will definitely be significant. 
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